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Abstract. We introduce a quasi-local integral functional and scalar quasi- 
local variables to examine a wide class of spherically symmetric inhomogeneous 
spacetimes that generalize the Lemaitre-Tolman-Bondi (LTB) dust solutions 
("LTB" spacetimes). By using these variables, we can transform the fluid 
flow evolution equations into evolution equations for non-linear, covariant, 
gauge-invariant perturbations of Priedman-Lemaitre-Robertson- Walker (FLRW) 
cosmologies. In the linear limit, we obtain spherical perturbations in the 
synchronous gauge under the long wavelength approximation. The formalism has 
a significant potential for cosmological applications, as it allows one to examine a 
wide variety of sources with different "equations of state", generalizing known 
FLRW solutions to idealized but non-trivial and non-linear inhomogeneous 
conditions. The quasi-local functional can be reformulated as a weighed proper 
volume average distribution, with the weight factor given by a scalar invariant 
related to the quasi-local mass-energy function. The back-reaction terms, 
emerging in Buchert's proper averaging formalism, can be expressed as differences 
between fluctuations of averaged and quasi-local energy densities. By comparing 
this average with the weighed quasi-local one, we can define a binding energy 
functional related to spatial gradients of the averaged and quasi-local variables 
that appear in the back-reaction terms. 



PACS numbers: 98.80.-k, 04.20.-q, 95. 36. -fx, 95.35.-fd 

1. Introduction. 

Recent observations seem to favor a dominant theoretical consensus, "the concordance 
model', in which present cosmic dynamics is dominated by an elusive source called 
"dark energy", somehow associated with a repulsive (yet unknown) interaction [1]. 
While the other dominant cosmic source, dark matter, is taken as inhomogeneous 
due to its association with structure formation at the galactic scale, dark energy is 
mostly considered as homogeneous, or nearly so, as it is assumed to have a signifhcant 
effect only at a larger cosmic scale in which the universe appears to be homogeneous. 
However, given our ignorance on the fundamental properties of these sources, there is 
no theoretically binding reason to assume, a priori^ that no valuable new information 
would result from studying their interaction under inhomogeneous conditions, at least 
in scales associated with structure formation and gravitational clustering. 

Since non-linear inhomogeneity in gravitational clustering or virialization of 
bound structure occurs at Newtonian velocities, these phenomena have been studied in 
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a Newtonian framework [5] |2l SI [S] . However, there is no reason to avoid using the low- 
energy regime of a relativistic theory, rather than a pure Newtonian framework. After 
all, rest-mass is the only source of Newtonian gravity, so it is unreasonable to assume 
that the main effect of these unknown sources can be reduced to rest-mass, ignoring 
possible effects from (yet unknown) forms of matter-energy contributions and their 
interactions. In fact, even with a cosmological constant, the study of the equilibrium 
and stability of astrophysical sources yields different results in the Newtonian limit of 
General Relativity from those of a pure Newtonian approach [6]. 

A large body of recent research considers a relativistic treatment of 
inhomogeneous sources in the galactic and cosmic scale as part of the critique of 
the concordance model. For example, the possibility that cosmic acceleration (or 
at least part of it) could be an effect of due to the presence of inhomogeneites in 
photon trajectories in observations from high red-shift objects [71 [8| l9l ITO l fTT | [T2 l [T3] . 
From a more theoretical angle, a number of articles relate or explain cosmic 
acceleration in connection with non-local effects that appear in various averaging 
formalisms [Ml [151 [TBI [TTl [181 [19] . In particular, an "effective" acceleration emerges 
in connection with the so-called "back-reaction" in Buchert's scalar average formalism 
[20] . This acceleration can be negative and could be, in principle, observed [21 ] [22 l [23] . 
Also, back-reaction can be associated with the existence of non-trivial gradients of 
quasi-local energy between gravitationally bound systems in an expanding cosmic 
background. [ni[2l[25| 

While the aim of most of the above mentioned proposals is to question the dark 
energy paradigm, they could also be compatible with it, as their predicted effects could 
provide just a partial explanation for the observed acceleration. However, even if the 
dark energy explanation prevails over the alternatives that dispute it, a controversy 
which is beyond the scope of this article, there is no harm done (and perhaps valuable 
new information) in probing its behavior under inhomogeneous conditions. As a 
compromise between "realistic" inhomogeneity that would require numerical methods 
and codes of high complexity and the approach of linear perturbations, we examine 
in this article inhomogeneous spacetimes under idealized but non-trivial conditions. 
We describe and summarize below the contents of the article. 

By assuming a spherically symmetric Lemaitre-Tolman-Bondi metric, we derive 
in section [2l a class of spacetimes ("LTB spacetimes") that generalize to nonzero 
pressure (and nonzero pressure gradients) the well know solutions for inhomogeneous 
dust associated with this metric [26l- The most general source for this metric (in 
the comoving frame) is a fluid with anisotropic stresses, which will end up being 
interpreted as fluctuations of the isotropic pressure. However, regardless of the 
interpretation, we remark that the existence of pressure anisotropy is far from a 
drawback or shortcoming, not only because it can be associated to a number of 
well motivated physical effects |27], but because an inhomogenous source with only 
isotropic pressure is a far more idealized situation than one with anisotropic pressure. 

In section [3] we show that LTB spacetimes can be fully characterized by a set 
of covariant scalars, whose evolution equations in the covariant fluid flow "1+3" 
representation [28l|29l[30] reduce to a system of non-linear autonomous scalar PDE's 
and constraints involving radial gradients. In section [4] we introduce a "quasi-local" 
integral functional, and its associated "quasi-local" scalar functions, related to the 
Misner-Sharp quasi-local mass-energy [3T1 [351 [331 [M] ■ Expressing all local covariant 
scalars as perturbations of the quasi-local scalars, we rewrite in section [5] the "1+3" 
fluid flow equations as evolution equations for the quasi-local variables and their 
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perturbations. These equations are effectively equivalent to a system of ODE's in 
which radial dependency enters as a parameter (see Appendix C of [55]). 

We prove in section [SI on the basis of known criteria that define a perturbation 
formalism [551 HE 133 IMj, that the quasi-local evolution equations can be 
rigorously characterized as evolution equations of covariant, gauge-invariant, non- 
linear perturbations on a FLRW formal "background" given by the quasi-local 
variables (which satisfy FLRW dynamical equations). All scalars associated with 
quantities that vanish in the FLRW background: Weyl tensor, anisotropic pressure 
and shear are expressible as gauge invariant first order fluctuations of the energy 
density, isotropic pressure and expansion scalar around their corresponding quasi- 
local variables. 

As we show in section [3 once we choose an "equation of state" (EOS) between 
the quasi-local energy density and pressure (which are "background" variables), the 
whole system of evolution equations becomes determined. Given this equation of 
state, we point out that non-local effects associated with the long range nature 
of the interactions in self-gravitating sources [37| can arise when comparing local 
and quasi-local scalars. In order to examine this issue, we look in section [8] at 
the classical ideal gas, showing that non-linear perturbations of particle number 
density and temperature provide a first order "virial" correction (given as a correlation 
energy) to the "usual" EOS satisfied by the background variables. This type of 
virial corrections to the ideal gas EOS have been found in Newtonian self-gravitating 
systems [SH [39l |40]. Similar effects could arise in applications of the perturbation 
formalism to self-gravitating relativistic dark energy sources whose fundamental 
physics is still unknown. 

We verify in section [S] that in the linear limit, when perturbations are "small", 
these evolution equations reduce to those of spherical linear perturbations in the 
synchronous gauge, but under the long wave length approximation [41] |42]. In 
section [10] we rewrite, in full covariant form, the quasi-local functional and scalar 
functions in terms of weighed proper volume integrals, where the "weight factor" 
is an invariant scalar associated with the quasi-local mass-energy. The quasi-local 
functional becomes then a weighed spatial average distribution that generalizes the 
standard spatial average. By comparing in section [Tl] the integral mass-energy 
functionals associated with the quasi-local (weighed) and standard (weight = 1) 
averages, we can define a domain-dependent binding energy. 

A theoretical connection between quasi-local variables and the issue of back- 
reaction has been suggested by Wiltshire [TT] |23J US] , in the context of the scalar 
averaging formalism of Buchert and coworkers. We provide useful arguments to this 
suggestion in section [T2l by showing that both kinematic and dynamical back-reaction 
terms are essentially comparisons of the square fluctuations of the expansion scalar 
and pressure around their spatial average and around their quasi-local equivalents. As 
shown elsewhere [43] , this comparison provides sufficient conditions for a positive back- 
reaction (and necessary for a negative "effective" acceleration) . These fluctuations can 
be expressed in terms of spatial gradients of the averaged and quasi-local variables, 
and also in terms of the domain dependent behavior of the binding energy. In section 
1131 we provide a brief summary of our results and suggest possible avenues for further 
research. We discuss several regularity issues in the Appendix. 

In follow up articles we will examine numeric applications of the formalism 
presented here, looking at the issue of back-reaction, as well as the modelling of 
cosmological dark matter and dark energy sources. 
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2. Generalized LTB spacetimes 

We will use the term "LTB spacetimes" to denote all solutions of Einstein's equations 
described by the spherically symmetric Lemaitre-Tolman-Bondi metric in a comoving 
frame: 

ds^ = ^C^de + + i?2[d02 ^ sin2 0^02]^ 

\ — K 

where R = R{ct,r), R' — dR/dr and K — K{r). LTB spacetimes comply with: 
(a) 4-velocity is a geodesic field, (b) zero energy flow in a comoving frame and (c) 
existence of a regular symmetry center [26^ . The metric ([ij is usually associated with 
a dust source, but the most general source compatible with it is an anisotropic fluid 

T'^^ = + (2) 

where u"- = (5Q,with ^{ct,r), p{ct,r) being, respectively, the energy density and 
isotropic pressure, H"'' is the anisotropic pressure trace-less tensor (H" — WaH"'' — 0) 
and hab ~ UaUjj+gab — ^a^ldij is the induced metric of hypersurfaces ^T(t) orthogonal 
to the 4-velocity {t constant, with space indexes i,j — r, 6, (p). 

Besides /x, p and IT''', the remaining covariant objects associated with LTB 
spacetimes are 

• Expansion scalar: 

• Ricci scalar of hypersurfaces ^T{t): 



3^ _ 2{KRy 



• Shear tensor: 



<^ab = V(qUq) - -hab, (5) 



• The "electric" Weyl tensor: 

E^'' = UelidC"^^^ (6) 

where R — dR/dct and (7°'"^'* is the Weyl curvature tensor. 

In LTB spacetimes all spacelike symmetric trace-less tensors like H"'', cr"'', E"-^ 
can be expressed in terms of a single scalar function as 

n-^fc^psah^ ^ ■P = S„^,^''^ (7a) 

„ab ^^^ab^ ^ I]==SahCT''\ (76) 

^ah^^^afc^ ^ £=^^^,E''\ (7c) 

where 5°'' = /i"'' — 3x"x'' and = \/hr^ S'jt is the unit vector orthogonal to w° and 
to the 2-spheres, orbits of SO (3). 

The scalars /x, p, V follow from the field equations G""^ — kT"'' (with k = SttG/c'*) 
for dl])-©: 

KpiR'^R' = R{R^ + K) ', (8a) 
KpR^R' = - i \R{Ff + K) + 2R^r] ' , (86) 
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while from ^ and (|7fl|) - ([7c|) . we obtain for £ and S 



E = l 
3 



R 
R 



R' 



R^ + K 



(9 a) 



(96) 



2 ■ 6 " ■ 2i?2 ■ 
It is evident that LTB spacetimes are completely characterized by the following set of 
covariant scalars 



{m, p, V, e, s, £, 37^}, 

which will be called the "local" scalar representation. 

The best known LTB spacetime is that with a dust source: p : 

pc^, while 



(10) 



= V = 0, for which 
reduce to 



R^ 



(11) 
(12) 



total energy density is rest-mass density: pi 

npc^R^R' ^ 2M', 
2M 

where M = M{r), which appears as an "integration constant', is the "effective" rest 
mass-energy in length units. The rest of the scalars 0, E, ^TZ, £ remain the same 
(with 7^ = in ((96|) ). The standard tactic to deal with the dust case is to solve the 
Friedman equation ([72]) for a given choice of M(r) and K(r), and then to evaluate 
from this solution all the derivatives of R, M, K in order to compute the scalars 
p, e, E, ^n, £ from ©, (gl), jMl and There is an extensive literature 

on how this is done [HI |45] (for an extensive review, see [26|), and the same tactic 
can work on simple particular cases with non-zero pressure and simple "equations of 
state" relating p, p and V (see [3S1I1S1I17]), but for a general LTB spacetime it is not 
possible (or very difficult) to find (not to mention solve) an equation like (I12|) in order 
to evaluate radial gradients of R for computing the scalars in ((TI 



3. The "1+3" decomposition 

An alternative approach to LTB spacetimes follows from the covariant "1+3" 
spacetime decomposition of Ehlers, Ellis, Dunsbury and van EUst [551 UHl 1301 j or "fiuid 
flow" equations, in which the field and conservation equations V^T"* = for ([I])-® 
are transformed into a first order system of time evolution equations and constraints 
for the covariant scalars pO| . 

Bearing in mind (|7ap -(l7c |) . the "1+3" equations for LTB spacetimes become the 
following set of scalar evolution equations 

e- -^-|(Ai + 3p)-6E2, (13a) 
p ^ ^ (/i + p)e-6E-p, (13&) 



26 



K 



E - -^E + E^-f + ^T', (13c) 
^ = ~^V~^[p + p-2V]J:-3[£ + ^v) (l + s), (13rf) 



together with the spacelike constraints 

IT 
R 



{p-2V)'^6V^^0, (14a) 
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E + l) +3E^=0, (14&) 

^(/^+^^)' + f' + 35f =0, (14c) 
and the Friedman equation (or "Hamiltonian" constraint) 

which is an integral of the Raychaudhuri equation (|13a|) . The system (|13ap - P^ 
requires an equation of state (EOS) Unking fj., p and V to become determined. 
However, time and radial derivatives do not decouple, in general. Moreover, pO| 
is not the only scalar representation for LTB spacetimes. 

4. The quasi local functional and scalar functions. 

The Misner-Sharp quasi-local mass-energy function is a well known invariant 
defined for spherically symmetric spacetimes as (HI [32l [33l |34] : 

—— = H K'- ia)ib)V(a)Vib)V V - J- : (-LOj 

^ 9rr 900 



where 'RS°-^'^^^ (c)(d) is the Riemann tensor in an orthonormal tetrad and i?^"^ = 
ip(b) — \/h^S^ are the tetrad unit vectors tangent to the orbits of S0(3). In fact, 
the two main invariants in spherically symmetric spacetimes are precisely M and R 
|33| , so that quantities expressible in terms of these invariants or their derivatives are 
necessarily covariant. 

For LTB spacetimes with an anisotropic source tensor like 1^ , the definition 
leads to: 

— =R^ + K, (17a) 

2M' = KfiR'^R\ (Ub) 
2M = - k{p-2V)R^R. (17c) 

whose integrability condition is equivalent to the conservation law VbT""^ = 
(equations (jl36p and (|14ap ). Notice that for the particular case of dust (^ = 7^ = 0) 
we have Ai = M and the quasi-local mass-energy is conserved along the 4-velocity 
flow: Ai = Q. When p, V are nonzero, then ()17ap generalizes the Friedman equation 

By looking at (|17aP and (Il76p , and assuming the existence of a regular symmetry 
center at r = 0, an integral expression for A4 can be constructed by integrating both 
sides of (|8al) along spherical comoving domains in hypersurfaces ^T{t): [33] 

^ i::oi^R'R'dx _ r^+k 

3 J^Z^R^R'dx i?2 • ^ ^ 

Different forms of mass-energy cannot be split or localized (in a covariant covariant) in 
General Relativity, but the function A4 obtained by means of (HH]) can be interpreted 
as a total or "effective" mass-energy associated with a compact domain bounded by 
the 2-sphere marked by r. This function also has a well understood behavior in the 
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appropriate asymptotic limits j3Il 1321 [Ml IMj ■ The integral ^TE\i motivates us to de- 
fine an integral functional acting on scalar functions defined in this type of compact 
domains: 



Definition 1. Quasi-local functional. Let X{T>) be the set of all smooth 
integrable scalar functions in the regular comoving domain I? = §^ x c 
^T{t), where §^ is the unit 2-sphere parametrized by {9, 0) and = {a; | < 
X < r}, with a; = marking a symmetry center 0. We define the "quasi-local 
functional" as the linear integral functional ( )* : X{T)) R, such that for 
every scalar function A G X{'D) we get the real number 

^ ' g R^R'dx ^ ' 

where the notation "[r]" emphasizes the fact that the dependence on r is 
that of a functional acting on functions in X{T>), to distinguish it from the 
dependence "(r)" in scalar function acting on T). In order to simplify nota- 
tion, we use Jg = J^^q to denote integration in the real interval and we 
will omit expressing explicitly the time dependence of A, R, R' and (A)*[r]. 
Unless specified otherwise, we will assume inside an integral symbol that: 
t is a constant parameter, all functions depend on x and R' — dR/dx. A 
covariant definition equivalent to will be given in sectionlTOl Meanwhile, 
we take this expression as an operative definition. 

Definition 2. Quasi-local scalars. We can use the functional ([19]) to construct 
scalar functions in T) with the same "correspondence rule" . Consider the map 
J, : X{V)~^ X{V), such that for every Ae X(V) and for aU ro € 

A, = J,{A):V^R, A,{ro) = {A)4ro]. (20) 
The real valued functions A^ will be denoted "quasi-local scalars" , as they 
generalize to any A G X the integral definition of the quasi-local mass- 
energy function (|18p constructed with fi. 

Basic properties. The scalar functions A* (r), comply with the following useful 
properties: 

1 f 

A'R^dx, (21a) 

R-^{r) Jo 

AJ ^^[A-A,], (21 &) 

-ft 

i, =(i)*+3(HA)*-3H,A, (21c) 
The functionals satisfy exactly the same properties above, with A* replaced 
by and taking [r] as an arbitrary parameter. In order to simplify 

notation, we have adopted the following conventions: 

A, EE [JM)]' ^ MA) = {A),, 

AJ = [JM)]' j^JM')^{^)*. 

(AB), = X{AB). 

(22) 

J Restrictions on T) occur if the ^T{t) have spherical topology or when there is a local curvature 
singularity. See the Appendix 
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It is straightforward to verify that these properties are perfectly self- 
consistent. 

Functional vs function inside integrals The quasi-local functional and scalars 
(|19p and ((20| behave differently under integration. The functions A* 
integrate as any scalar function, hence if we apply ( ), to A* we get: 

C AJx)R^R'dx 

whereas, (yl)»[r] is just a real number that depends on the domain of a 
previous integral (|19p . hence it complies with the following iteration rule: 

which allows us to regard the functional (A),[r] (but not the function A^,{r)) 
as an integral average distribution on A. See section [TT] 

Basic quasi-local scalars. Considering (fT8|) and the definitions (fT9|) and (|20|) . 
the quasi-local matter-energy energy /i* = J7* (/i) is related to the quasi-local 
mass-energy function by 

3^^* = ^. (25) 

From ([3]), (|4]) and evaluating p9)) - ([20l) for the scalars 9, '^TZ, their quasi-local 
equivalents of are 

e* =J,(e) = ^, (26a) 
H 

In the following sections, until section [TOl we will deal only with quasi-local 
scalar functions 

By applying ([Ill)-(I2ni) to (Eill-lHg), we get 

K J nR^R'dx = ^R{k^ + , (27a) 

K j pR^R'Ax = '^P*R^ " "3 ^ (-^^ + -ft^) + 2i?i? , (276) 

which, with the help from (|26ap - (|266l) . lead to equations for p,, 0, and ^7?., that 
are identical to the Raychaudhuri and Friedman equation of FLRW models: 

e, = -^-|(/i*+3p,), (28) 

= 3/.*-—, (29) 

where 8j — [J7*(8)]^ 7^ J7*(0^). These equations can be combined to yield the same 
FLRW energy conservation equation 

/i* = — (^i* + P*) B*. (30) 



Quasi-local variables, non-linear perturbations and back-reaction 



Applying (fTn)) - PD|) to (|5c|l. (Pa|l and (|9&| yields the remaining scalars, T,,V and £, as 
deviations or fluctuations of O, /i and p from their quasi-local equivalents 



E = --[e-e,], 
^ ^ [p~p*], 



£ = - 



6 



(31a) 
(31&) 
(31c) 



while the time derivative of the quasi-local mass-energy complies with 

2M = -Kp^R^R. (32) 

Since A4 and R are invariants of spherically symmetric spacetimes, it is evident that 
/It*, and 0* are covariant scalars, related to these invariants and their derivatives 
by means of (l25|), ((26^ and p2l) . 



5. Evolution equations for the quasi local scalars. 

A very convenient relation between the scalars A and A^, follows by defining the 
"relative deviations" or "perturbations" 
A A^ 



5iA) 



A = A^ 



1 + <5(-4) 



(33) 



which leads to the alternative covariant scalar representation of LTB spacetimes given 

by 

{e.,/.,,p./7^.,<5(«),<5(^),<5(p),<5(^^)}. (34) 

We will call (|34p the "quasi-local" scalar representation and remark that it is fully 
equivalent to the local representation (fTU)) . Regularity conditions for the 5^"^^ functions 
are discussed in the Appendix. 

The "1+3" system (I13ap - P^ is the set of evolution and constraint equations 
for (jlOp. we derive now the evolution and constraint equations for the representation 
(|34p. From (j21bl) . the radial gradients of /i,, and can be given in terms of the 
S functions by 

while (pS)) and ([50)1 are evolution equations for /i* and 8*. Hence, the evolution 
equations for S'^^'^ and (J^^^ follow from the consistency condition 



[AJ]-. 



(36) 

applied to ([28|l . ([30]) and ([35]) for A = 0* and /i*. The result is the following set of 
autonomous evolution equations for the quasi-local scalar representation (134^ : 

[1 + w] 6,, (37a) 



= - 



a, 

3 



(l + 5(«)) 



5(e) 



6(e,/3) 



(37&) 
(37c) 
(37rf) 
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where 




(38) 



The constraints associated with these evolution equations are simply the spatial 
gradients ([35|l . while the Friedman equation (or Hamiltonian constraint) is ((29|) . 
Notice that the constraints ([55)1 follow directly from differentiating the integral 
definition (fT9)) . so by using the quasi-local variables we do not need to solve these 
constraints in order to integrate p7ap - (|37(ip . 

With the help of (|31ap - (|31cp and (|33p . it is straightforwards to prove that 
the evolution equations (|37ap - (j37(ij) and the constraints (|29p and ([55)1 are wholly 
equivalent to the 1+3 evolution equations (|13ap - (|13(j|) and their constraints (|14ap - 
(|14cl) and (fTSj) for the local scalar representation pO|) . Hence, they are also equivalent 
to the field and conservation equations. The system (|37ap - (j37dl) becomes fully 
determined once an "equation of state" that fixes as a function of /i* is selected 
(see section [7]) • 

6. A non linear perturbation scheme 

The dimensionless functions J^"^' in ([55)1 allow us to express the local scalars A as 
non-linear relative fluctuations or perturbations of their quasi-local equivalents A^. 
Notice, from (l^^ - (|266p and (HH!)-®, that the scalars {6,, ^7^*} satisfy 

FLRW evolution laws. However, these are just the mappings under (|^D|) of the LTB 
covariant scalars {Q, ^i,p, ^TZ}, whose FLRW equivalents are also covariant scalars. 
This suggests that the role of the (5^"^' as spherical perturbations on a formal FLRW 
"background" state given by the A^ can be rigorously justified by means of the known 
gauge invariant and/or covariant perturbation formalisms developed by Bardeen [36] 
and Elhs, Bruni and Dunsbury [IHl HSl [50] . 

A perturbation formalism (not necessarily linear) between an idealized spacetime, 
S (a FLRW cosmology), and an "lumpy" universe model, S, follows by defining a 
"background model" in S, constructed by objects (in 5") that are the images of a 
suitable map $ between objects in S to objects in S |2H]. Though, perturbations in 
general are not uniquely defined by such an abstract map <&, because of the ambiguity 
of the gauge freedom associated with choices of coordinates and hypersurfaces. This 
fact requires carefully defining perturbed variables that are "gauge-invariant" [SBj. 
However, in the specific case of setting a perturbation formalism when S is an LTB 
spacetime we face a simpler task, as both the idealized FLRW and the perturbed LTB 
spacetimes can be completely described by scalar functions. Hence, $ can be a map of 
scalars to scalars, while the fact that both spacetimes are spherically symmetric and 
are given in the same normal geodesic coordinate (or frame) representation, greatly 
suppresses most of the gauge freedom that we would find for a general S. 

For the FLRW- LTB case under consideration, the map <i> can be defined rigorously 
as follows. Let X and X be, respectively, the sets of smooth integrable scalar functions 
in S and S, then for all covariant FLRW scalars A ^ X (we denote FLRW objects 
with an over-bar) the map 



defines a "background model" (associated to a FLRW cosmology) in LTB spacetimes 
through the quasi-local scalars A^ (which are LTB objects satisfying FLRW 



$ : A -> A, 



<i>{A)^X{A)^A,eX, 



(39) 
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dynamics). Given their common normal geodesic representation and considering the 
perturbations 

S^^^ - (40) 

associated with (pQ)) . LTB spacetimes in the metric ([T]) with source ([2]) can be 
considered then as spherical non-linear "perturbed" FLRW cosmologies in the 
"synchronous" gauge. 

The perturbation scheme that we described above can also be understood in 
terms of the gauge invariant and covariant (GIG) approach of Dunsbury, Ellis and 
Bruni [551 [23]. Following these authors, a perturbation scheme on FLRW cosmologies 
is covariant if the "lumpy" model S is described by variables defined in the framework 
of the "1+3" fluid flow variables /i, p, 8 of (|13ap - (|13dl) . Although our description of 
LTB spacetimes is not based on these scalars, it is still a covariant description because 
/X*, p* and are covariant scalars by virtue of their connection with the invariants 
A4, R and their derivatives in ([25]) . (|26a|) and ([32| . Hence, the formalism is covariant. 

As commented by Ellis and Bruni [28], by virtue of the Stewart-Walker gauge 
invariance lemma [48] . all covariant objects in S that would vanish in the background 
S (a FLRW cosmology in this case) are gauge invariant (GI), to all orders, and also 
in the usual sense (as in [36]). The basic tensorial quantities in LTB spacetimes that 
vanish for a FLRW cosmology are the tensors H"'', cr°^ and E'^'' given by ([2]), ([5]) and 
([6]), which from ([7a|) -([7c j) and (|31a|) - picp . can be written as 

^ [9 - e*] (41a) 



3 



1 



6 

1 

ab 



3 , 



'^ab 



(Alb) 



U^" = -[p-p.]E'^'', (41c) 

where is the tensor defined in (|7ap -(|7c p . which is (for spherical symmetry) a 
covariant object. From (|216|) . ([55)1 and ([iTal) - (|irc|) . it is evident that 

— H*, p — p* and — 0* and the radial gradients p'^ and Q'^ are all "first 
order" quantities (in the perturbation scheme) that are GI to all orders. Hence, from 
(|33|) . the perturbation variables S^^\ S^p^^ and S'^^^ are also GI to all orders, though 
/L(*, p*, 0* are not (which is expected because these are "zero order" background 
variables). Hence, the fluid flow dynamics of LTB spacetimes in the quasi-local scalar 
representation ([M)) is that of spherical, non-linear GIC perturbations on a FLRW 
background. In section [9] we examine the connection with linear perturbations. 

It is worth remarking that ^TZ^, does not appear in the dynamic equations ([37aj)- 
(|37rf|l - though it can appear when specifying initial conditions to integrate these 
equations. This variable is GI only if the FLRW cosmology S is spatially flat, but its 
associated variables '^'^ and ^TV^ are GI. In fact, from the Friedman equation ([25)1 
and pS)) we obtain 



6 =3^*-lTj ' ^^2") 
!|i5(^^) = ^;.,^(M)-2(|l)%(^), (426) 
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so that we can always eliminate either one of 8* or 5'--^\ (S'®-* in terms of '^7?.* and 
and construct a system of evolution equations equivalent to p7a|) - (|37(ij) . but 
describing the dynamics in terms of spatial curvature ^TZ^ and its perturbation 5^ 

7. Equations of state. 

The system p7ap - (j37dp is still undetermined, as there are no evolution equations for 
and S^P\ Equations (|37ap and (j376p are formally identical to FLRW equations, 
and so are the evolution equations for the background, while (|37cp and (j37dp are the 
evolution equations for the perturbations and thus convey the effects of inhomogeneity. 
As in any perturbative approach, we need to impose an "equation of state" (EOS) 
between and /i, to determine the background subsystem (j37ap - (j37bp . Such a choice 
of an EOS also determines the perturbation equations. 

Consider a commonly used "barotropic" equation of the form 

= (43) 

From (j35p . we obtain the corresponding relation between fluctuations and 
perturbations of /i and p 

dp* , 1 ,,, . 

P-P* = j^il^- l^-*h (44a) 

^(^) ^dlnp,^(^)^ ^^^^^ 
cfln^* 

An important and simple particular case of (|43p is given by the linear EOS 

= Wo M* 1 (45) 

where Wo — Wq (t) , and in particular wq can be a constant (so that = is dust) . In 
this case we have 

p-p,=wo[^i-^i,], 5(p)=5(a'). (46) 

Notice that, in general, if an EOS like ([43|l is imposed it will not be satisfied by 
p and fi. But this is to be expected, since in a perturbation scheme fluctuations or 
perturbations do not satisfy, in general, the EOS of the background (not even in linear 
perturbations). The exception is the linear EOS ([43]), which implies p — wq ^, even if 
Wq = wo{t) because the quasi-local functional (fT9)) involves integrals in a hypersurface 
of constant t, hence J*{wQ{t) fi) = wo{t) J,(^) for all t. 

It could be argued that, even if a given is a physically reasonable EOS, it 
could yield a wholly unphysical relation between p and /x through (|44ap . While this 
possibility cannot be ruled out, it is something that must be tested by looking at 
the solutions of (|37al) - (j37dj) for the particular EOS. If such a situation arises, then 
the specific LTB spacetime associated with this EOS could fail to provide an adequate 
non-linear perturbative description for the FLRW cosmology with this particular EOS. 
However, even if such a development could occur for a given EOS, there is no reason 
to expect in self-gravitating systems that an EOS between /x* and p* should also hold 
between local variables p and /i. This expectation is only correct in homogeneous 
conditions or for hydrodynamical systems characterized by short range interactions, 
much smaller than macroscopic scales [37] ■ Therefore, there is no reason to demand 
or to force such behavior of the EOS if we apply our formalism to self-gravitating 
relativistic sources of cosmological interest, like dark energy or the Chaplygin gas. 
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whose fundamental interactions are still unknown and cannot be assumed to be of a 
hydrodynamical nature. 

Evidently, the plausibility of the non-linear perturbation formalism must be 
tested and judged for different EOS according to its predictions. In this context, 
it is worth remarking that quasi-local variables convey non-local information because 
they are integral quantities that depend on the behavior of the local variables in the 
whole integration range D = x ^r- Since gravitation is a long range interaction, this 
non-local information can be relevant in the dynamical description of self-gravitating 
systems. We illustrate this fact in the following section using as example the classical 
ideal gas. Also, it has been suggested [171 [Ml [25] that "back-reaction" is related to 
dynamical and observational effects of quasi-local energy. We discuss this issue in 
section |l2l 



8. Quasi local variables in the classical ideal gas. 



In order to illustrate directly the perturbation formalism under a physically motivated 
and familiar EOS, we consider the classical ideal gas. This example illustrates the type 
of effects that could arise if we apply this scheme to sources like dark matter or dark 
energy. Following the arguments of the previous section, instead of (|43p we choose the 
classical ideal gas EOS for background variables /i* and 



(47) 



where — Ji,{n) and T* = J*{T) are, respectively, the quasi-local scalar functions 
associated by means of (|20p to particle number density and temperature, and k is 
Boltzmann's constant. The EOS (|47l) leads to a relation between the perturbations 
and (5("\ 5'^^\ so that n, T relate to n*, by means of The choice 

of background EOS (|47l) transforms (j37a|) - (|37rf|) into a set of evolution equations 
for n*, T*, 5^"^ 5^'^\ 6* and 6'^^\ leading to a description of this source as an 
inhomogeneous self-gravitating system in the non-linear perturbation formalism 
under consideration. 

It is important to mention that the same EOS (|T7)) . but relating local quantities 
n and T, is strictly valid for non-interacting particles, or as an approximation to 
hydrodynamical sources under short range interactions (see |46j for a hydrodynamical 
approach to the ideal gas by means of LTB spacetimes). For a self-gravitating system, 
given in terms of non-linear perturbations from an ideal gas in a FLRW background, 
we do not expect that fi and p will be related by (|47l) . 

From (|35|). ([33]) and ([47]) we obtain the perturbations of /i and p in terms of 
perturbations of n and T 



S^^\ (48a) 

(486) 

Using ([551) (I48al) - (|486p we can express fi and p in terms of the local variables n, T 
and their fluctuations as 



mc n 



-nkT 



— mc^n H — nkT 



(n — n*) (T — r*) 

<5(")5(T) 



(H-<5("))(l + (5(^)) 



(49 a) 
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P 



nkT 



1 



(n 



(T - n) 



uT 



nkT 



1 



(l + (5("))(l + ^(^)) 



(496) 



which shows that (|T7|) leads to forms for the local pressure and internal energy 
[(3/2)nfcT] that are equivalent to the "standard" (or "non-perturbed") EOS plus 
a virial correction that depends on the correlation between fluctuations of n and T. 
This can also be appreciated by using the integral form (|21ap . leading to 

{n-n^){T-T^) _ f n'R^dx J T'R^dx 



nT 



nTR^ 

J J dxdx R^{ct, x)R^(ct, i) (dn/dx) (dT/dx) 



(50) 



nTR^ 

which is clearly a correlation density analogous to those found for Newtonian self- 
gravitational systems ^38j [39l [40] . Although these articles considered gravitationally 
bound systems, and so their correlation terms will be different from those in a free- 
falling LTB ideal gas model, the qualitative comparison that we make with them 
is relevant, as it shows the effect of prescribing an EOS like (|T7)) on the relation 
between "perturbed" variables, fi andp, under a non-linear perturbation scheme. This 
perturbative treatment furnishes non-local information that adds correction terms to 
the standard ideal gas EOS "p — nkT", which is strictly valid for non-interacting 
particles. This extra information can be related to the long range interaction effects of 
gravity (and to "back-reaction", see section [T^ in scales where (n — n*)(T — r*)/(nT) 
is not negligible. Since fluctuations like {n — n^){T — T^) are related to the gradients 
of and T^, by (j216p . the virial correction terms in (|49a[) - (j496p are important when 
there is significant spatial variation in n and T. 



9. Connection with linear perturbations. 

Given the fact that S'-^^ 5'-p\ (5^^) describe LTB spacetimes as spherical non-linear 
perturbations with respect to the "background" variables /z*, p*, O*, it is important 
to examine the evolution equations p7ap - (|37dj) in the linear regime, when these 
fluctuations are "small", i.e. S'^^^ ^1. In these conditions we should be able to 
treat them (at least formally) as linear perturbations. 

In order to make this comparison, we derive a second order equation for the 
density perturbation (5*^^^ for the special case of a barotropic EOS (|45| with wq 
constant. By differentiating both sides of (|37aP and using the remaining equations 
in (|37ap - (j37rf|) to eliminate all derivatives except (5*^^^ and S'^^\ we obtain 

'^^'^ - , , + ^ ^^'^ - S (1 + 3"'o)/i* S(^^ h + wo + -5(^)1 = 0, (51) 

1 + Wo + o*-^' 3 2 L J 

which is an exact equation. For near homogeneous conditions, as assumed in linear 
perturbations, ([STjl reduces to 

SM + ^1 sM - I (1 + 3wo){l + wo)^i. 5^''^ = 0, (52) 

an equation that is formally identical to the evolution equation for linear density 
perturbations of a fluid with this EOS around a FLRW background (/x*, 6*) in the 
synchronous gauge |4T], and under the conditions of perturbation amplitudes much 
larger than the horizon [HI Notice that in the case of dust, wq = 0, we recover 
the well known perturbation equations in the comoving gauge (which for dust is 
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a synchronous gauge as well) and without any restriction on the amplitude of the 
perturbations. 

It is straightforward to show (with more complicated algebraic manipulations) 
that, for all EOS among quasi-local variables, the evolution of 6'-^'^ also reduces in the 
linear limit to linear density perturbation equations in the synchronous gauge and in 
the long wavelength limit. This is not surprising since the metric ([T]) of LTB spacetimes 
is given in coordinates that are synchronous (because of its geodesic 4-velocity) and 
comoving. There are known technical difficulties with an equation like (|52p in this 
gauge (see Appendix of ^?^). However, the inconsistencies shown in this reference do 
not apply to the linear limit of LTB spacetimes with w ^ becuse the latter are also 
comoving and so — 0" from the outset in equations (A1)-(A3) of [42]. Thus, the 
linear perturbation limit in the synchronous gauge with long wavelengths is consistent. 
In a non-linear regime the correct equation is (j51|l . not ([52]) . therefore the problems 
with this gauge in connection with the latter equation do not apply, as S^"^^ <C 1 is no 
longer valid. 



10. A covariant quasi— local functional. 

A covariant expression for the quasi-local functional (fT9|l and scalar functions ([20]) 
follows rewriting these definitions in terms of integrals over the proper spatial volume 
element dV ~ y^det^hab) d'^x, where d^x = dr d^ d0. Considering the same domain 
P = X C ^T{t), with £,r — {x\0 < X < r}, the proper spatial volume for V in 
LTB spacetimes is 

V = / dV = 47r / R^R'dx, (53a) 

T=[1-KY'^, (53&) 

where the Att comes from the integration over §^ . The definitions (|19p and become 
"weighed" proper volume integrals: 

where 

dV* = TdV, V* = / TdV = ^R^. (55) 

The "weight factor" (|53b[) is a conserved invariant quantity (i.e. = 0), which by 
virtue of (|17a[) . can be given as 



i?2 



R J 



(56) 



Notice that T is an invariant because R = u°'VaR, and M and R are invariants in 
spherically symmetric spacetimes [HI [32l |33l [34] . As mentioned in section [H A4 can be 
interpreted as a total "effective" mass-energy in V. In the flat space limit Ai/R ^ 0, 
this invariant reduces to the Lorentz contraction factor [49]. Also, if we replace A4 
in (|56p with the "Schwarzschild mass" this equation is identical to the equation for 
radial geodesies in Schwarzschild geometry, with J-^ playing the role of the conserved 
quantity associated with the timelike Killing field (related to the binding energy). 
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We show in the foUowing section that J- is related to a binding energy integraL 
Another intuitive interpretation for follows from the Newtonian limit [33]. The 
leading term of A4 in this limit is the Newtonian rest mass integral: Ai « 
(G/c2 )Af„ewt, while R = {l/c)dR/dt < 1 and « GM„ewt/(c^i?) < 1, hence 

(l56l) becomes at first order 



1 



dR 
'dt 



GM„ 



R 



(57) 



which is = 1" plus a correction of order given by the sum of Newtonian kinetic 
and gravitational potential energies per mass unit for comoving shells. In a relativistic 
context this splitting of energy contributions is no longer valid. 



11. Quasi local vs proper averages. 

The quasi-local functional, as expressed in covariant form by ([5^ . can be considered 
a weighed average distribution with the conserved invariant playing the role of a 
formal "probability density" function. By its construction as a functional and its 
domain dependence, ([54]) complies with the properties of average distributions of 
random variables, such as the iteration rule: ((A),), = {A)^ given in ([Ml) , as well 
as the definition of second order momenta (variance and covariance): 

{{A-{A),r), ^{A'),-(A)l (58a) 

{{A - {A),){B - (B),), = (AB), - {A)4B),, {58b) 

where we have omitted the [r], as the functional domain dependence is clear. The 
quasi-local averages constructed with (jM]) satisfy the mathematical properties piap - 
(pT^ : 

A-{A), = -i^ / A'V.dx, (59a) 
V*{r) Jo 

(A) J (596) 

(A),- = {A),+3{nA), -3(n)4A),, (59c) 

where (A),' = do{A)^, ^ {doA)^, = (v4)*. It is evident from ([23]) that quasi-local scalar 
functions A^, do not satisfy (|58ap - (|58&[) . and so they are not "averages" of any sort. 

The quasi-local average is a generalization of the standard or proper ("weight" 
= 1) spatial average functional: ( ) : X{'D) — > R, such that for every A e XlV) we 
get the real number 

LAdV AT-^R^R'dx , ^ 

The average functionals ( ) and ( ), satisfy the same properties of average distributions 
of random variables, (|24)) and l|58ap - (|586p . Also, the proper average satisfies (|59ap - 
(|59cl) . with (A)* and V, replaced by (A) and V. The form of the radial gradients of 
(A) in (|59fe[) expressed in terms of metric functions is: 

(^)' = V - (^)] = - (^)] ' (61) 
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where we applied the foUowing relations: 



y =(-^>, (62a) 
(A). = {62b) 



that follow directly from ^ and ([60)1 . 

Considering M as the "effective" mass-energy in V given by the integral of in 
(|176[) . as discussed in scctionSl suggests applying (f54|) and (l60|) to /x. This leads to the 
following quasi-local and proper mass-energy functionals, related with the averages 
of fi and their corresponding volumes by 

2M[r] = / /idV, = ^(/i)*V*, (63a) 
Jv c 

2Mp,[r]= / fxdV^^{fi)V. (636) 
Jv c 

(63c) 

As shown by these two expressions, in general, quasi-local ("effective") and proper 
mass-energy is different and this difference (in spherical symmetry) can be used to 
define a "binding energy" |50| . which in our case is the functional 

fiW-— [(m>*V,-(a.)V] = — / [T-l]^,AV^ — {{T-l)^,)V, (64) 

C C J-p c 

where we used (I62a[) - (|62fep . We have then an interpretation for JF, as an invariant 
scalar function such that, for a given energy density /i and domain V, the sign of 
— I determines the sign of its corresponding binding energy B functional (|64p . 
Since pressure is also an energy density, similar mass-energy and binding energy 
functionals can be associated with p 
2G f 2G 

W*W = ^ / pdV, = ^(p)*V*, (65a) 
2G f 2G 

WprM = — / pdV= — (p)V, (656) 
c- Jv c 

on r Of 
[r] = - Wpr - — / [T~l]pdV^—{{T-l)p)V, (65c) 
c Jv c 

Clearly, and Wpr can be interpreted as the proper and quasi-local work done by 

the averaged pressures (p)* and {p) in the domain T). In fact, comparing ((65a|) with 

([32]), we get: 2M = ->V*V* = -W*V, (6),. 

Regarding the averages (8)* and (O), they are simply the proper time rate of 

change of V* and V: 

V, it, V a 

where we have introduced the "scale factor" a oc V^/'^[r] to characterize the (domain 

1/3 

dependent) length scale of the averaging ( ), just as i? (x V* [r] (curvature radius or 
"area-distance" in V) characterizes the domain dependent scale of ( Thus, for a 
dust source R{r) oc {p)V^[r\ and a(r) cx {p)^/^[r\. 

From Ip ^ -(|626 p - (p5)) . it is clear that the ratios 



(m) Mpr ip) Wpr 
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depend on the invariant J-', and so, on the binding energy functional and on the spatial 
curvature through the function K — {^TZ^ R^)/6 in (|26&p . However, this is a non-local 
dependence, as it is determined by the behavior of the functions /i and K in V. If 
we have positive, zero or negative spatial curvature in the whole 2?, then J- — I and 
B will be (respectively) negative, zero or positive. However, if a domain D contains a 
local region with strong positive curvature (undergoing local collapse) , binding energy 
can be still be positive if '^TZ is asymptotically negative and D sufficiently extended. 
Hence, as opposed to an FLRW cosmology, in an inhomogeneous LTB model the 
binding energy is not only domain dependent, but can change signs from one domain 
V to another in each '^T. If spatial curvature vanishes in T), then = 1 and ( = ( ), 
binding energy is zero and the quasi-local averages become just proper averages. 

As shown in [171 EQI Ell ES] , the scale dependence of binding energies (associated 
with A4 and A^pr) are related to the so-called "back-reaction", and could provide 
important information connected with the theoretical interpretation of current 
cosmological observations. We examine this issue in the following section. 



12. Connection with "Back reaction". 



"Back-reaction" terms appear in a scalar averaging formalism, such as that developed 
by Buchert [50], when the proper volume average functional (pO)) is applied to local 
scalar functions in fluid flow evolution equations like (|13a[) - (jl3(ip . As we show in 
this section, the same terms arise in LTB spacetimes with both functionals (|54[) and 
(|60| . While for LTB spacetimes the fluid flow evolution equations, as well as p7a|) - 
(j37(j|) . are scalar covariant equations, for general spacetimes, the only covariant scalar 
equations are the Raychaudhuri, energy conservation and Friedman equations, (|13ap . 
(jl36|) and We rewrite these equations for LTB spacetimes as 

e = -^-|[M + 3p]-^(e-e.)^ (esa) 
/i = -(M+p)e + (e-e,)(p-p,), (68&) 

where we used piap - (|31cp to express S and V in terms of fluctuations of and p 
around their quasi-local equivalents 8* and p*. 

Considering the derivative commutation property (|59cP and (|58ap - (j586p . we apply 
the spatial average functional (|60p to (|68ap - (|68cp and work out the terms so that the 
averages (0), {p) satisfy evolution equations that are FLRW plus "corrections". 
The result is 

(6)- = -M!-^[(^)+3(p)] + Q«, (69a) 
(m)- = -[(A^) + (p)](e)-Q('*\ (696) 

where the "corrections" Q^'^^ and Q'^^^ are, respectively, the kinematic and dynamical 
"back-reaction" terms, given by 

Q('^) = ^(cM), cM = (e-(e))'-(e-e,)', (70) 
sw = (c(^)), c('i)^(e-(e))(p-(p))-(e-e.)(p-p.), (7i) 
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which coincide with the terms found elsewhere [20j . As a comparison, if we apply 
the quasi-local average functional to (|68ap - (|68cp . we obtain evolution equations 
identical to (|69a|) - (|69cp . with ( )^ replacing ( ), and with back-reaction terms given 

by 

s?^ = ^(cW)*, d"^) = (e - (9), )' - (e - e,)' , (72) 
Ql^^ = (d^^)„ d''^ = (e-(e>,)(p~(p),)-(e-e,)(p-P*), (73) 



which are exact analogues of ([70)1 and ([7T|l above. Because of the fact that A»(r) = 
{A)4r] (from the terms ci^^ and ci'^^ in ^ and ^ vanish at X — r, but this 

would not make qI"^^ = and qI*^^ = because for cc < r, both C^'' and C»'^'' are 
nonzero, and so the average ( )* would integrate over them (over x) yielding nonzero 
back-reaction terms. Only if spatial curvature is zero in all of T) we would have 
= = in all P, and then there would be no back-reaction [43]. If using 
( ), back-reaction also vanishes when spatial curvature is zero in all 2?, since then 
(> = (>*■ 

The basic motivation for using back-reaction terms in Buchert's formalism is to 
obtain an "effective" negative acceleration, Acs_, that would mimic cosmic acceleration 
detected by observations. Considering a dust source: p = (p) =P* = Q*'*) = and 
jjL = pc^ , where p is rest-mass density, the condition A^s < follows by rewriting 
(|69al) - (|69cp in terms of A^s, leading to the following condition (see ^20J for details) 

Aeff - f (P) - Q*'^ < 0, (74) 



where the kinematic back-reaction term, Q^'^-', is given by ([70)1 . 

As shown in [l^, a sufficient condition for Q'"^^ > (and necessary for ([74]) ) is 
simply 

C^^) > 0, (75) 



in the whole domain V. Clearly, ([75]) indicates that the necessary condition for ([741) 
is basically a comparison between fluctuations of 6 with respect to its average, (9), 
with fluctuations with respect to its quasi-local equivalent scalar, 9* (but not the 
quasi-local average (9)*). In LTB spacetimes with nonzero pressure, we must add 
the "dynamical" back-reaction with the comparison of fluctuations p — (p) and p — p^, 
in ([ZH). 

It has been suggested by Wiltshire [TTl dS] that back-reaction, and so cosmic 
acceleration, is related to the spatial gradients of quasi-local energies between bounded 
structures that have decoupled from an expanding cosmic background. It is evident 
that this should be the case, since the fluctuations 9 — 9, and 9 — (9) that appear 
in ([70]) . ([TT]) and ([75[) are related to radial gradients of 9, and (9) by means of ([35)1 
and ([H]) 

9', =^[9 -9,] = ^[9 -9,], (76a) 

which translates into gradients of quasi-local energy, since (because of (126 ap ) 9, — 
3R/R can be identified with a quasi-local "kinetic" energy density in the left hand side 
of Friedman equation ([29)) . in its form ([56)) . and in the Newtonian approximation ([57]). 
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while (Q) can be understood as an average "kinetic" energy density in this context. 
Since J- and {J-') are directly related to the binding energy in (|64p . we can relate 
these gradients to the bound states of observers in different domains V. In fact, the 
back-reaction terms can be directly written as averages of gradients. Using psp and 
(|6ip . we rewrite (iTOt and ((7T|) as a comparison of spatial gradients: 



whose value, by virtue of (|62a|) - (j62&|) and ((64|) . depends on how JF determines the 
change of sign of the binding energy in the averaging domain V. These results are in 
agreement with and can give (model dependent) support to the arguments given by 
Wiltshire in [171 HH HI] • The connection between quasi-local variables and functionals 
with the issue of back-reaction deserves a proper examination that is being undertaken 
in a separate article. 

13. Conclusion and possible avenues for future research. 

We have examined a class of spherically symmetric inhomogeneous spacetimes 
(denoted "LTB spacetimes") associated with the LTB metric ([1]), whose source is 
an anisotropic fluid ([2]). In general, LTB spacetimes can be compatible with a wide 
variety of sources and "equations of state" , and in particular, we can always recover 
the well known LTB dust solutions by setting p — 0. 

Following the theoretical formalism that we have presented in this article (which 
we summarized in the introduction), LTB spacetimes can be used as models of 
interacting dark matter and dark energy, as well as models to test alternative proposals 
based on back-reaction effects. As non-linear perturbations of a FLRW cosmology 
in the quasi-local scalar representation, they are ideal theoretical tools to construct 
appealing models in cosmology that generalize (to idealized but non-trivial and non- 
linear inhomogeneity conditions) a large ammount of known solutions in modern 
cosmology that have been derived and examined in a FLRW context (or with linear 
perturbations). 

While the integration of the evolution equations (j37a|) - (j37rfj) definitely requires 
numerical work, these equations can be handled as a system of ODE's in which radial 
dependence enters as a parameter (see Appendix C of [35 thus the evolution system 
can be treated as a "dynamical system" using numerical techniques that are quite 
accessible and can be handled by known and widely used packages. In fact, the quasi- 
local scalar representation has already been used successfully to conduct a "dynamical 
system" study of LTB dust solutions. [35] 

The connection with back-reaction is straightforward. The "kinematic" and 
"dynamical" back-reaction terms that emerge from Buchert's averaging formalism 
are both expressible as differences between squared fluctuations (variances and 
covariances) of & and p around their proper averages, (8) and (p), and their quasi- 
local equivalents, 0* and p*. These fluctuations can also be given in terms of 
spatial gradients of (9), (p) and 6*, p*, and these expressions can also be related 
to the (domain dependent) binding energy integral B. These comparisons provide a 
quantitative (though model dependent) framework to derive the conditions for the 




(77) 



(78) 
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existence of back-reaction [13], and a subsequent negative "effective" acceleration 
that could mimic cosmic acceleration. This framework provides also a substantative 
appraisal to the the qualitative arguments of Wiltshire [T71 [Ml that characterize 
back-reaction as an effect related to the existence of spatial variation in quasi-local 
energy between bound observers and an expanding cosmic background. 

Evidently, the predictions of LTB spacetimes, as sources that are inhomogeneous 
at specific scales, must be tested in reference to observations [II [3 [H [TOl [HJ [T2j [13] 
(see also [5T1[23] and references quoted therein). However, at least in scales below the 
so-called homogeneity scale (100-300 Mpc), models built with LTB spacetimes can 
provide useful information about non-linear effects associated with inhomogeneity, 
such as "spherical collapse" models depicting gravitational collapse in local regions 
embedded in an expanding homogenous background, effects of dark energy in large 
scale dynamics of galactic superclusters, back-reaction terms in the dynamics of 
averaged scalars, quasi-local energies and non-local phenomena. In a follow up article 
we will apply LTB spacetimes as models for inhomogeneous dark energy sources, such 
as the Chaplygin gas, interactive mixtures of dark matter and dark energy and scalar 
fields. We will also re-examine in more depth, and in separate articles, the connection 
between quasi-local variables and back-reaction (and the "effective" acceleration). 

It is important to stress that the quasi-local representation can be reformulated 
in terms of fluctuations, variances and covariances of quasi-local average distributions 
that can be compared to those associated with the "standard" proper volume average 
of Buchert's formalism. In particular, the quasi-local average introduces a "weight 
function" that is an invariant that could be interpreted as some sort of non- 
trivial probability density. Further developments of his framework could lead to a 
definition of a quasi-local "gravitational entropy" functional (in the context of [52] ). 
and/or to possible theoretical analogies with statistical mechanics and black hole 
thermodynamics. This is a very appealing possible avenue for future research and 
is presently under consideration. 

Finally, the formalism that we have presented can be readily generalized beyond 
spherical symmetry to the "quasi-spherical" Szekeres models, which do not admite 
isometries |26[ 155] . While the metric coefficients, as well as all covariant quantities, 
in these models depend on all coordinates, it has been proven that the proper spatial 
volume average of scalars is independent of the "non-radial" coordinates [S^ (see also 
[55]). As we prove in a forthcomming article, this is also the case with the quasi-local 
functions and averages. Thus, the techinques used in this article are valid for these 
general models with just minimal adjustments, but the latter models are no longer 
spherically symmetric and so provide a less idealized framework. 

Appendix: some regularity issues 

As mentioned in section IV, we have assumed that the integration domain I> = x 
with = {x \ < X < r}, containing a symmetry center at x = 0, is fully regular. In 
particular, if the hypersurfaces '^T{t) have spherical S'^ topology, there is a second sym- 
metry center, x — r^, and so the domain is restricted to a; < re- Another restriction 
comes from the well known fact that curvature singularities (initial and collapsing) 
arise in LTB dust models. In particular, there is a re-collapsing singularity in regions 
with positive spatial curvature, K > 0, and when 9 < [311I3S]. Since this type of 
curvature singularities is bound to occur as well in the more general LTB models, it 
is worthwhile commenting on its possible effects. 
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The coordinate locii of initial or collapsing singularities in dust solutions are (in 
general) not simultaneous {i.e. not marked by a constant t). Even if parameters can 
be set up to have a simultaneous initial singularity [Ml US], in general, the collaps- 
ing singularity in regions with K > is not simultaneous, but marked by a curve 
[ct{r coll), rcoii] in the (ct,r) coordinate plane, so that R{ct{r coii),rcoii) = 0, and cur- 
vature scalars diverge as comoving observers approach these coordinates (notice that 
this is different from a symmetry center, r = rc, where R{ct,rc) — regularly for all 
t 7^ t{rcoii)). Thus, in any such collapsing region the hypersurfaces ^T{t) for t > t(rcoii) 
are only regular for the semi open subset ^coii = {x \ TcoIi < x < r} C ^r- However, 
the existence of this (and the initial) singularity has no consequence in the definition 
of the quasi-local functional and scalar functions, because the involved integrals can 
be treated simply as standard improper integrals. We define at each ^T(t) the incum- 
bent integrals with their lower integration limit asy = TcoIi -t-e, for an arbitrarily small 
e > 0, and then obtain the limit as e ^ 0. Off course, since Q — > — oo in this limit, 0, 
and (6), might diverge as well, but these functions are well defined in the range ^coU- 

Another important regularity issue, also present in dust LTB solutions, is the 
possibility of unphysical "shell crossing" singularities, which might occur if i?' = 
for non-comoving coordinate values {i.e. coordinate values not equal to constant 
r). See [44l |45] for a comprehensive discussion. It is important to distinguish the 
shell crossings from "turning values" at a given constant value of r = rtv where 
R'{rtv) = regularly, as in the case where the ^T{t) have spherical topology (or 
for "closed" FLRW spatial sections). The possible effect of a shell crossing singularity 
in LTB spacetimes would be to make the functions J*^"^^ ill defined. From (j21b[) . ([55)1 
and (|596p . the S^"^^ will be regular if i?' = implies A'^ = 0, which will happen at 
r = rtv But if R'{ct,r) = and r ^ rtv, then R' = will happen with A'^ ^ 0, and 
consequently, 8^'^'> diverges. Since A = A,[l-|-(5(^)], this would mean that A diverges 
for a finite A*, a wholly unphysical and unacceptable situation. Therefore, a standard 
regularity condition in LTB dust solutions is to demand absence of shell crossings, or 
R! > 0, save at turning values for hypersurfaces with spherical (S'^) or "wormhole" 
topologies (S^ X R or X §^) that have no symmetry centers. This regularity condition 
should also be demanded for the more general LTB spacetimes, which could exhibit 
spherical topology. For the "wormhole" topologies, the integrals ([T9[) and ([54[) defining 
the quasi-local functional could be unbounded, but this could be avoided if the lower 
integration limit is defined asymptotically (see [35j for an example). 
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